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PREFACE 


Personnel  of  the  Environmental  Systems  Division  (ESD) , Mobility 
and  Environmental  Systems  Laboratory  (MESL),  U.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES),  conducted  this  study  for  the  U.  S.  Army 
Materiel  Development  and  Readiness  Command  in  furtherance  of  DA  Project 
No.  1T162112A528,  "Environmental  Effects  on  Materiel,"  Task  02,  "Terres- 
trial Effects  on  Materiel  Development,"  and  in  support  of  the  Project 
Manager,  Remotely  Monitored  Battlefield  Sensor  Systems,  under  Support 
Agreement  No.  AMC-PM-F.BS- 74-06-18. 

The  study  was  conducted  during  the  period  August  1973  - September 
1974  under  the  general  supervision  of  Messrs.  W.  G.  Shockley,  Chief, 

MESL,  and  B.  0.  Benn , former  Chief,  Environmental  Research  Branch  (ERB) 
and  now  Chief,  ESD,  and  under  the  direct  supervision  of  Mr.  J.  R.  Lundien, 
former  Chief,  ERB,  and  now  Senior  Electrical  Engineer,  ESD  Project 
leader  was  Dr.  D.  H.  Cress.  The  computer  programs  used  in  the  report  were 
compiled  by  Messrs.  J.  R.  Stabler  and  E.  A.  Baylot,  ERB.  The  report  was 
prepared  by  Dr.  Cress.  Discussions  with  Messrs.  Carl  Bennett  and  Raymond 
Fidler  of  the  Naval  Coastal  Systems  Laboratory,  Panama  City,  Florida, 
concerning  analogous  location  techniques  using  pressure  waves  in  water  were 
very  helpful. 

Director  of  WES  during  this  study  and  preparation  of  the  report  was 
COL  G.  H.  Hilt,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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V.  S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


Multiply 

millimetres 
centimetres 
metres 
kilometres 
metres  per  second 
kilometres  per  hour 


: used  in  this  report 

can  be  converted  as  follows 

By  . . . 

To  Obtain 

Metric  (SI)  to  U.  S. 

Customary 

0.03937007 

inches 

0.03937007 

inches 

3.280839 

feet 

0,6213711 

miles  (U.  S.  statute) 

3.280839 

feet  per  second 

0.6213711 

miles  (U.  S.  statute)  per 

U.  S.  Customary  to  Metric  (SI) 

tons  (short) 
inches  per  second 
degrees  (angular) 


907.1847 

25.4 

0.01745329 


kilograms 

millimetres  per  second 
radians 


4 


SEISMIC  METHODS  OF  LOCATING  MILITARY  GROUND  TARGETS 


PART  I : INTRODUCTION 


background 


1.  In  its  effort  to  acquire  the  capability  to  detect  and  locate 
ground-based  targets,  the  Army  has  developed  two  concepts,  one  in  which 
an  array  of  unattended  seismic  sensors  is  deployed  in  a region  of  sur- 
veillance to  detect  and  locate  continuously  moving  targets  (men,  vehi- 
cles, etc.)  and  the  other  in  which  acoustic  sensors  are  placed  around 
the  region  of  surveillance  to  locate  firing  artillery.  Moving  targets 
are  located  by  monitoring  the  activations  of  the  individual  sensors  in 
the  array,  the  accuracy  depending  upon  the  range  at  which  the  sensors 
are  activated.  The  individual  sensors  are  spaced  such  that  the  activa- 
tions allow  quasi-continuous  tracking  of  the  target  without  an  over- 
lapping of  the  ranges  of  activation  (truly  continuous  tracking  by  this 
method  is  not  possible).  Thus,  a large  number  of  sensors  are  necessary. 
Targets  generating  impulsive  signatures,  such  as  firing  artillery,  car. 
be  located  using  three  or  more  sensors.  The  time  of  arrival  of  the 
acoustic  signature  at  each  sensor  is  noted,  and  the  artillery  is  located 
by  uniquely  relating  the  differences  in  these  arrival  times  to  the 
coordinates  of  the  target.  The  difference  in  arrival  times  between  any 
two  sensors  can  be  used  to  construct  hyperbolic  lines  of  possible  target 
positions.  Intersection  points  for  the  constructed  lines  of  several 
pairs  of  sensors  mark  the  most  likely  position  of  the  target.  A problem 
arises  if  the  ranges  from  multiple  firing  locations  to  the  sensors  are 
such  that  the  ordering  of  the  arrival  of  the  acoustic  signature  at  a 
given  sensor  is  different  from  the  ordering  of  the  arrivals  at  other 
sensors.  The  sensitivity  of  the  instrumentation  to  the  intensity  of  the 
signals  and  the  sensitivity  of  the  algorithms  to  errors  in  the  arrival 
times  are  such  that  the  region  of  surveillance  must  be  surrounded  by 
sensors.  Thus,  current  target-location  systems  require  that  the 
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target  be  within  the  area  bounded  by  the  arrays  of  sensors  (seismic  or 
acoustic)  and  that  the  region  of  surveillance  be  accessible  to  per- 
sonnel installing  them. 

Approach  and  Scope  of  This  Study 

2.  In  a search  for  improved  methods  of  using  seismic  energy  to 
locate  ground  targets,  this  study  has  been  directed  toward  evaluating 
the  feasibility  of  using  a target-location  system  consisting  of  a group 
of  geophones,  i.e.  three  geophones  configured  such  that  a geophone  is 
placed  at  each  of  the  vertices  of  an  equilateral  triangle  whose  sides 
are  oriented  to  a specified  reference.  The  geophones  are  spaced  less 
than  7 m apart;  this  array  will  be  referred  to  herein  as  a short-based 
array.  A target  can  be  located  in  a two-dimensional  space,  as  shown  in 
Figure  1,  provided  the  range  from  some  reference  point  (sensor  location) 
to  the  target  and  the  directional  angle  (i.e.  the  angle  between  the  line 
from  the  sensor  location  to  the  target  and  some  reference  direction)  are 
V wu.  There  are  several  concepts  that  may  be  used  for  determining 
target  location,  as  will  be  discussed  later.  However,  the  directional 
angle  from  one  or  more  short-based  arrays  must  always  be  obtained.  By 
determining  the  phase  shift  in  the  seismic  wave  train  received  at  each 
geophone,  the  direction  to  the  target  can  be  calculated.  Since  seismic 
waves  in  the  ground  and  acoustic  energy  that  may  couple  with  the  ground 
at  the  geophone  travel  at  much  higher  velocities  than  do  most  land 
vehicles  (targets) , the  travel  distances  for  moving  targets  will  be 
small  during  the  time  the  signal  travels  from  the  target  to  the  geo- 
phone. Therefore,  if  processing  times  are  kept  short,  the  same  tech- 
niques used  for  a stationary  target  c^rv  be  used  for  a moving  target. 

3.  The  use  of  short-based  arrays  of  geophones  for  determining 
direction  from  an  array  to  a target  has  not  previously  been  rigorously 
investigated;  however,  the  concept  has  been  used  successfully  for 
tracking  targets  that  generate  pressure  waves  in  water.  Furthermore, 
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Figure  . Parameters  (i.e.  directional  angle  and  range) 
required  for  location  of  a target  in  a two-dimensional 

space 


the  Naval  Coastal  Systems  Laboratory*  has  demonstrated  that  an  array  of 
hydrophones  can  be  used  to  track  several  different  energy  sources 
simultaneously,  thus  offering  the  hope  that  the  approach  can  be  used  to 
locate  multiple  targets  using  signatures  in  seismic  media  as  well. 

4.  The  measure  of  feasibility,  as  treated  in  this  report,  is  the 
accuracy  attainable  with  short-based  arrays  of  geophones  in  locating 
targets.  Other  factors  affecting  the  development  of  the  system,  such  as 
cost  and  deployment  requirements,  are  not  considered.  However,  the 
level  of  sophistication  being  pursued  for  other  location  systems  (for 
example,  the  bigradient  microphone  technique  at  the  Naval  Ordnance 
Laboratory)  is  equivalent  in  complexity  to  the  signal  processing  and 
instrumentation  specifications  necessary  for  implementation  of  the 
techniques  reported  herein. 

5.  The  accuracy  with  which  the  direction  from  the  array  to  the 
particular  target  can  be  determined  depends  on  the  degree  to  which  the 
seismic  wave  fronts  propagate  in  a straight  line  from  the  target  to  the 
geophones,  the  response  characteristics  of  the  instrumentation,  and  the 
accuracy  with  which  the  signature  from  the  particular  target  can  be 
separated  from  the  signatures  arising  from  background  noise  and  other 
targets.  Therefore,  the  following  factors  have  been  considered  within 
the  scope  of  this  study: 

ci.  Signal-processing  requirements  for  implementing  the 

target-location  technique  in  single-  and  multiple-target 
environments. 

_b.  The  relation  between  design  parameters  for  the  short- 

based  array  (i.e.  spacing  between  geophones,  spatial  and 
phase  reponse  of  geophones,  and  phase  response  of 
amplifiers)  and  accuracy. 

6.  To  address  the  factors  listed  above,  two  concepts  of  target 
location  are  presented  in  Part  II,  both  of  which  require  determination 
of  the  directional  angle.  A directional  algorithm  is  derived  for  the 
short-based  array.  Two  techniques  of  signal  processing  for  application 


* C.  M.  Bennett,  ''The  Directional  Analysis  of  Ocean  Waves,"  Informal 
Report  NCSL  144-72,  Dec  1972,  Coastal  Systems  Laboratory,  Panama  City, 
Fla. 
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of  the  directional  algorithm  are  presented;  they  differ  in  their  degrees 
of  complexity.  In  Part  III  these  techniques  are  applied  to  test  data 
obtained  with  wheeled  and  tracked  vehicles  in  two  test  programs  at  two 
geologically  dissimilar  sites  (Vicksburg,  Mississippi,  and  Fort  Carson, 
Colorado).  Conclusions  and  recommendations  are  presented  in  Part  IV. 

A thorough  evaluation  of  the  sensitivity  of  the  directional  algorithm 
to  terrain  variations,  such  as  terrain-induced  changes  in  the  propaga- 
tion path  of  the  seismic  energy,  was  beyond  the  scope  of  this  study. 
However,  other  sources  of  error,  spatial  and  temporal,  are  discussed 
in  Appendix  A.  Appendix  B presents  a possible  technique  for  locating 
targets  by  considering  the  differences  in  arrival  times  of  seismic  and 
acoustic  waves.  Vehicle  test  conditions  are  tabulated  in  Appendix  C. 


PART  II:  CONCEPTS  FOR  TARGET  LOCATION 


7.  A target  can  be  located  by  obtaining  the  range  and  directional 
angle  relative  to  a reference  point  and  reference  direction  (Figure  1) . 
The  reference  point  can  be  defined  as  the  location  of  a short-based 
array  of  seismic  sensors.  The  target  can  then  be  located  by  several 
schemes.  Among  the  concepts  of  location  are: 

a_.  Location  by  triangulation  in  which  the  directional  angles 
to  a particular  target  from  two  spatially  separated 
short-based  arrays  are  determined  (Figure  2).  The  point 
of  intersection  of  the  two  lines  defined  by  the  direc- 
tional angles  is  the  estimated  location  of  the  target. 
Triangulation  uniquely  determines  the  two  coordinates 
required  for  location  because  the  directional  angles  are 
determined  from  the  data  received  at  the  reference  array, 
and  the  range  is  uniquely  defined  by  the  distance  from  the 
reference  array  to  the  intersection  of  the  lines  along  the 
directional  angles  determined  by  the  spatially  separated 
arrays . 

b^.  Ranging  from  a short-based  array  to  a target  by  monitoring 
the  difference  in  time  of  arrival  of  two  modes  of  energy 
transport  (such  as  seismic  and  acout tic  energies)  having 
differences  in  their  respective  velocities  of  propagation. 
This  technique  of  target  location  appears  applicable  only 
to  impulsive  targets,  such  as  firing  artillery  and  impacts. 
It  is  discussed  farther  in  Appendix  B. 

It  should  be  noted  that  the  results  of  the  study  reported  herein  are  more 
applicable  to  target  location  by  triangulat j.on  than  to  target  location 
by  ranging,  since  the  data  analyzed  in  this  report  consist  of  signatures 
of  moving  targets  rather  than  impulsive  (short-lifetime)  signatures. 
However,  the  ranging  technique  is  included  in  Appendix  B because  it  ap- 
pears to  he  a possible  approach  to  location  of  impulsive  sources  based 
upon  a preliminary  examination  of  seismic  and  acoustic  data. 

8.  Both  of  the  approaches  to  target  location  described  in  the 
previous  paragraph  require  that  the  directional  angles  relative  to  two 
spatially  separated  short-based  arrays  (paragraph  7a)  or  one  short-based 
array  (paragraph  7b)  be  determined.  The  feasibility  of  using  short-based 
arrays  to  determine  the  location  of  a target  can,  to  a certain  degree,  be 
reduced  to  evaluating  the  accuracy  with  which  a short-based  array  can  be 
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used  to  determine  the  directional  angle  to  the  target.  Therefore,  the 
primary  effort  of  this  study  was  such  an  evaluation. 

9.  The  accuracy  with  which  a short-based  array  can  be  used  to 
obtain  the  directional  angle  can  be  described  by  some  angular  range  of 
indeterminacy  (A0)  (Figure  3).  The  A0  is  probabilistic  in  nature  and 
could  arbitrarily  be  defined  as  the  angular  range  for  which  the  total 
probability  of  the  target  being  within  this  range  is  less  than  a 
standard  deviation  from  the  actual  location.  For  purposes  of  analysis, 
it  is  net  necessary  to  statistically  define  the  A6  . The  region  of 
indeterminacy  (AR)  may  be  partially  or  totally  dependent  upon  the 
uncertainty  in  the  directional  angle  calculation.  If  the  target  is 
located  by  triangulation,  as  illustrated  in  Figure  3,  the  AR  is 
totally  defined  by  the  A 6 . Therefore,  the  accuracy  of  location  by 
triangulation  can  be  treated  by  examining  the  errors  in  the  directional 
angle  calculations.  However,  the  AR  resulting  from  target  location 
with  a single  short-based  array  and  ranging  by  time  of  arrival  differ- 
ences (paragraph  7b)  depends  upon  both  the  A9  and  errors  in  calcu- 
lation of  the  range.  An  approach  for  effectively  factoring  out  terrain- 
induced  errors  in  ranging  and  in  the  A0  for  stationary  or  quasi- 
stationary targets  is  discussed  in  Appendix  B;  a rigid  analysis  of  the 
errors  ip  ranging  is  beyond  the  scope  of  this  study. 

Derivation  of  Directional  Algorithm 

10.  As  the  surface  of  a medium  (soil,  rock,  etc.)  is  disturbed  by 
a forcing  function,  such  as  a vehicle,  an  artillery  piece,  or  a man, 
energy  is  propagated  outward  from  the  point  of  the  impulse.  The  mode  of 
propagation  can  be  Rayleigh  (i.e.  surface)  waves,  transverse  waves,  or 
corapressional  waves.  The  extent  to  which  cylindrical  symmetry  (i.e. 
radial  symmetry  about  an  axis  normal  to  the  surface  through  the  poi  t of 
the  impulse)  is  retained  by  the  wave  front  depends  upon  seismic  uni- 
formity of  the  propagating  medium.  The  medium  is  seismically  uniform  if 
the  velocity  of  propagation  is  not  spatially  dependent.  The  amount  of 
displacement  of  the  medium  from  equilibrium  is  a function  of  the  range 
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Region  of 


LEGEND 

Region  of  indeterminacy  (AR) 
Short-based  array 
Figure  3.  Definition  of  region  of  indeterminacy 
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from  the  impulse  to  the  location  at  which  the  medium  is  observed.  The 
displacement  can  be  in  the  radial  direction,  or  in  either  of  the  two 
transverse  directions  (i.e.  displacements  perpendicular  to  the  direction 
of  propagation) . The  amount  of  displacement  from  equilibrium,  for 
radial  or  transverse  motion,  is  denoted  by  £(r,t)  , where  r is  the 
range  to  the  impulse  and  t is  the  time.  Since  the  propagation  is 
assumed  to  be  radially  outward,  the  amount  of  displacement  arising  from 
a single  frequency  source  can  be  written 

S(r,t)  - A(r)  exp  j j <o  £t  - + 9 

where 

A(r)  = the  amplitude  of  the  wave 

j - *cr 

a)  * angular  frequency 
v(w)  « velocity 

6 = arbitrary  phase  angle 
The  variables  in  parentheses,  r and  u>  , denote  e functional  dependence 
on  the  respective  variables. 

11.  In  the  following  derivation  and  application  of  the  directional 
algorithm,  the  range  from  the  geophones  to  the  source  of  the  energy  is 
assumed  to  be  sufficiently  large  so  that  the  waves  can  be  approximated 
as  plane  waves  (i.e.,  the  two-dimensional  wave  front  defined  by  the 
spatial  distribution  of  the  wave  having  constant  phase  is  spatially 
described  by  straight  lines  rather  than  by  its  actual  circular  shape). 

As  shown  in  Appendix  A,  the  important  quantity  for  describing  the 
accuracy  of  these  assumptions  is  the  ratio  of  the  range  r to  the 
geophone  spacing  Z in  the  short-based  array.  Within  the  Remotely 
Monitored  Battlefield  Sensor  System  (REMBASS)  minimum  accuracy  require- 
ments, a ratio  of  11  or  greater  is  acceptable.  Within  the  area  con- 
tained inside  the  short-based  array  for  which  the  plane  wave  approxi- 
mation is  valid,  the  amplitude  of  the  wave  can  be  assumed  to  be  indepen- 
dent of  the  range.  A rigorous  treatment  of  the  possible  error  arising 
from  the  plane  wave  assumption  is  presented  in  Appendix  A. 
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12.  The  plane  waves  approach  the  short-based  array  (geophones  A, 
B,  and  C)  as  shown  in  Figure  4.  The  quantity  sensed  by  the  geophones  is 
the  particle  velocity  5 , which  is  obtained  by  taking  the  time  deriv- 
ative  of  £ . The  resulting  expression  for  S is  given  by 


The  assumption  of  plane  waves  means  that  the  wave  fronts  defined  by  a 
constant*'  phase  form  a line  perpendicular  to  the  direction  of  propaga- 
tion. In  Figure  4,  the  range  from  the  source  to  geophone  A is  denoted 
by  r^  . Since  the  wave  fronts  are  perpendicular  to  the  direction  of 
propagation  through  r^  , the  important  quantity  for  defining  the  phase 
difference  between  the  waves  at  two  points  (for  example,  at  geophones  A 
and  B)  is  the  projection  of  the  distance  between  the  two  points  upon 
the  direction  of  propagation.  This  projection  for  geophones  A and  B is 

JL_  cos  (3  - where  is  the  distance  between  A and  B,  (£  - 

AB  AB 

is  the  angle  between  the  direction  of  propagation  and  the  line  from  A to 
B,  and  </>  is  the  directional  angle.  Consistent  with  the  approximation 
of  plane  waves,  the  range  from  the  source  to  geophone  B is  given  by 

rB  ‘ rA  “ *AB  C0S  (P  " (2) 

Similarly,  the  range  from  the  source  to  C is  given  by 

rC  “ rA  ■ *AC  C0S  * (3) 

Equations  2 and  3 are  consistent  with  the  plane  wave  approximation 
because  their  accuracy  rests  upon  the  approximation  that  the  directions 
of  propagation  of  the  waves  passing  through  all  three  geophones  are  the 
same.  This  is  true  only  if  the  waves  are  accurately  described  as  plane 
waves  across  the  space  of  the  array. 

13.  Equations  2 and  3 can  be  substituted  into  Equation  1 to 
obtain  the  expressions  for  particle  velocity  at  geophones  A,  B,  and  C, 
as  follows: 
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SA  » wA  exp  j 


1 {“  [c  - vt] 


+ ® + f 


= wA  exp  j 


Sc  « ooA  exp  j 


'{“[*- 

{■[“ 


rA  “ *AB  COS  (g  ' ** 


vCw) 


■]+  9 + i) 


r A - & . „ C03  <|>  “1 

v(w)  ]*  9 + f 


Both  Equation  4 and  Equation  5 can  be  written  in  the  form 


^k=B,C 


r r 

1“ 

t - A 

v(m) 

+ 0 + 1 + A6Ak 


(4) 


(5) 


where  A 6^  is  the  phase  shift  of  the  signatures  at  B and  C relative  to 
A.  Therefore,  the  phase  shifts  are  given  by 


A0 


AB 


-u)iAg  cos  (3  - <{>) 

v(o>) 


(6) 


and 


-oj^Ac  cos  $ 


A0 


AC 


v(w) 


(7) 


The  time  increment  between  the  instant  B receives  the  signal  and  the 
instant  A receives  the  signal  is  given  by 


AB 


£ab  cos  (3  - <f>) 
v(u>) 


(8) 


Similarly,  the  time  increment  between  the  instant  A receives  the 
signal  and  the  instant  C receives  the  signal  is  given  by 

= *AC  COS  1)1 
TAC  v (a)) 

Cos  (3  - <J>)  can  be  expanded  using  a trigonometric  identity.  The  expanded 
form  of  Equation  8 can  then  be  divided  by  Equation  9.  The  resulting 
expression  for  <j>  is  found  to  be 


(9) 
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for  0 < 4>  < 360°  CIO) 


(j>  = 


-1 


tan 


r i /^actab 

Lsin  13  V*abt*c 


COS 


14.  The  spatial  symmetry  is  optimized  if  the  three  geophones  are 
placed  at  the  vertices  of  an  equilateral  triangle.  3 is  then  60  deg 
and  Equation  10  reduces  to 


4>  = 


for  0 < 41  < 360° 


(11) 


Note  that  the  signs  of  x^g  and  x^,  must  be  used  to  determine  the 
proper  quadrant  for  <j>  . 

15.  In  the  derivation  of  Equations  10  and  11,  a single  frequency 
source  was  assumed.  However,  the  forcing  function  of  the  source  is, 
in  general,  impulsive  in  nature.  Therefore,  a wide  band  of  frequencies 
is  introduced  into  the  medium.  The  particle  velocity  at  the  geophone 
is  described  by  the  continuous  sum  of  all  frequencies  according  to  the 
expression 


£(r , t) 


CO 

/• 


ojA(io) 


exp  j 


+ 0 


+ 1 
2 


doo 


The  signature  can  be  processed,  frequency  by  frequency,  for  directional 
estimates  by  obtaining  x^g  and  x^c  at  each  frequency,  x^g  and  x^, 
are  functions  of  frequency  because  the  wave  velocity  is  frequency 
dependent.  However,  as  long  as  the  frequency  increments  chosen  for 
directional  analysis  are  sufficiently  narrow,  such  that  the  wave  velocity 
in  a given  frequency  increment  can  be  considered  to  be  constant,  applica- 
tion of  Equations  10  and  11  to  obtain  directional  estimates  is  justified. 


Methods  of  Signature  Analysis 

16.  The  signatures  received  by  the  geophones  in  the  short-based 
array  must  be  processed  in  some  manner  to  obtain  the  time  delays  x^g 
and  x^  (Equations  8 and  9).  Equivalently,  the  phase  shifts  A0^g 
and  A0^£  (Equations  6 and  7)  can  be  obtained  for  each  frequency  w . 


18 


As  stated  in  paragraph  15,  the  signatures  generally  contain  a broad 
spectrum  of  frequencies.  The  degree  of  difficulty  in  analyzing  the 
signatures  is  roughly  proportional  to  the  width  of  the  spectral  band 
that  is  processed  to  obtain  and  . For  seismic  signatures  of 

military  targets,  the  spectral  region  of  interest  is  generally  between 
0 and  400  Hz.  As  the  width  of  the  spectral  band  being  processed  is 
narrowed,  the  capabilities  of  the  signal-processing  technique  are 
enhanced.  This  fact  becomes  particularly  apparent  in  multiple-target,  or 
noisy-background  environments  because  the  "target"  are  likely  to  trans- 
mit power  differently  for  different  spectral  regions  of  the  signature 
as  indicated  in  Figure  5a  for  two  targets  (1  and  2) . Since  the  signatures 
of  targets  1 and  2 are  superimposed,  the  power  spectrum  observed  at  the 
geophone  array  has  the  form  indicated  in  Figure  5b.  The  superposition 
of  the  signatures  of  targets  1 and  2 presents  a problem  in  interpreting 
the  calculated  time  increments  and  tac  . In  the  shaded  frequency 

region  of  the  spectrum  (Figure  5b) , the  signatures  of  both  targets  make 
approximately  the  same  contributions  to  the  spectrum  and,  therefore, 
have  approximately  equivalent  "rights”  to  the  resulting  calculated  value 
of  tAB  or  tAC  * Generally,  the  calculated  directional  angle  will 
fall  somewhere  between  the  two  targets.  However,  in  the  spectral  regions 
where  the  power  spectrum  of  target  1 dominates  over  that  of  target  2, 
the  calculated  values  of  and  result  in  accurate  calculations 

of  the  direction  to  target  1.  In  a practical  situation,  the  observer 
has  no  way  of  knowing  from  the  calculated  directional  angles  alone 
which  angles  are  accurate  descriptors  of  the  direction  to  the  target (s). 
Other  properties  of  the  signatures  that  describe  the  degree  of  super- 
position of  the  waves  must  be  invoked  to  aid  in  this  decision,  which 
adds  an  additional  complexity  to  the  signal-processing  technique. 

17.  Although  several  techniques  are  available  for  signal  processing, 
two  techniques  have  been  selected  that  have  differing  degrees  of  com- 
plexity and  corresponding  differences  in  their  capabilities.  The  two 
techniques  will  be  referred  to  as  cross-correlation  analysis  and 
coherent— signature  analysis.  The  cross-correlation  analysis  was  con- 
sidered because  of  its  simplicity.  However,  because  of  the  shortcomings 


19 


Power 


a.  Power  spectra  for  signatures  emanating  from 

two  targets 


b.  Power  spectrum  of  signature  at  the  array  in  a 
two-target  (targets  1 and  2)  environment 


Figure  5.  Contributions  of  signatures  of  two  targets  to 
the  power  spectrum  observed  at  the  array  in  a two-target 

environment 
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of  cross-correlation  analysis  (discussed  later) , the  coherent-signature 
analysis  was  implemented.  The  coherent-signature  analysis  has  a demon- 
strated applicability  to  processing  of  pressure  waves  in  water  for  cal- 
culating directional  angles  using  an  array  of  pressure  transducers.* 

The  primary  difference  between  the  two  analysis  techniques  (i.e.  cross- 
correlation and  coherent-signature)  is  the  frequency  bands  in  which 

and  (Equations  8 and  9)  are  computed.  The  cross-ccrrelation 


AB 


analysis  yields  values  of  and  xAr,  in  the  frequency  region  of 


AC 


maximum  energy;  the  coherent-signature  analysis  results  in  calculated 
values  for  and  x in  a number  of  narrow  frequency  bands. 

18.  It  should  be  noted  that  the  requirement  for  continuous  signa- 
tures over  time  for  coherent-signature  analysis  excludes  its  applica- 
tion, as  implemented  in  this  study,  to  impulsive  targets  such  as 
artillery  firing;  however,  cross-correlation  analysis  is  applicable  to 
impulsive  signatures. 

Cross-correlation  analysis 

19.  Cross-correlation  analysis  is  based  on  the  cross-covariance 

estimate  **  C _ for  a time  shift  u between  two  signatures  and 


where 


T/2 


cij(u)  * Y q(t)  + u) 


dt 


(12) 


-1/2 


The  cross-covariance  is  estimated  for  each  pair  of  geophones  i and  j 

in  the  triangular  array  and  for  a series  of  time  shifts  u such  that  u 

is  less  than  T . T is  the  length  of  the  integration  interval  and  must 

be  greater  than  both  the  period  of  the  dominant  frequency  band  and  the 

anticipated  time  delay  between  signatures,  and  dt  is  an  increment  in 

the  time  domain.  The  times  t and  t._,  are  then  obtained  by  finding 

AiJ  AC 

the  values  of  u for  which  C._  and  C1/n  are  maximized.  It  should  be 

AB  AC 

emphasized  that  the  times  and  x^,  are  frequency  dependent  because 


* Ibid.,  p.  8. 

**  G.  M.  Jenkins  and  D.  G.  Watts,  Spectral  Analysis  and  Its  Applications, 
Holden-Day,  San  Francisco,  Calif.,  1969. 
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of  dispersion,  and  the  times  obtained  using  cross-correlation  analysis 
are  characteristic  of  the  frequencies  of  the  dominant  energy  band. 
Coherent-signature  analysis 

20.  In  principle,  it  should  be  possible  to  compute  directional 
angles  for  each  frequency  present  in  the  seismic  signature.  In  prac- 
tice, however,  direct  analysis  using  a discrete  Fourier  transform  of  the 
raw  data  results  in  a scattering  of  the  computed  values  of  x and 
taC  and  a resulting  scattering  of  the  computed  directional  angles  for 
each  frequency.  One  approach  to  diminishing  the  amount  of  scattering  is 
to  average  the  x and  x^c  or  parameters  related  to  their  computa- 
tion (i.e.  Fourier  coefficients,  etc.)  over  time.  Although  an  averaging 
technique  is  incorporated  into  the  coherent-signature  analysis,  it  does 
not  sufficiently  reduce  the  scattering  in  the  calculated  values  of 
and  xAC  , primarily  because  the  length  of  the  averaging  time  cannot  be 
longer  than  several  seconds  for  meaningful  directional  calculations  to 
moving  targets.  Therefore,  some  criterion  is  needed  that  will  separate 
the  potentially  accurate  values  of  and  tac  from  the  remaining 

ones.  The  property  of  the  signatures  that  describes  the  consistency  in 
the  calculated  phase  shifts  between  the  geophones  in  the  array  has 
proven  useful  in  evaluating  the  accuracy  of  the  computed  value  of  x ^ 
and  x^c  . The  parameter  associated  with  this  property  of  the  signa- 
tures is  the  coherence.*  The  coherence  estimate,  as  determined  for 
the  signature  analysis  reported  subsequently,  has  the  following  properties: 

a_.  If  the  power  of  the  signature  remains  reasonably  constant 
at  a given  frequency  during  the  signature  time  required 
for  estimating  coherence  (approximately  2 sec),  and  the 
phase  shift  remains  constant,  the  coherence  is  approxi- 
mately unity. 

Ik  If  the  phase  shift  changes  with  time,  the  value  of  the 
coherence  becomes  progressively  smaller  with  increases  in 
the  magnitude  of  the  time-dependent  changes  in  phase 
shift  and  the  leng  *>f  the  averaging  time. 

It  seems  reasonable  to  assume  tha^  i..  r continuous  seismic  signatures 

emanating  from  a target,  the  phase  shift  should  be  reasonably  consistent 

over  relatively  short  time  intervals  (approximately  2 sec)  and,  therefore, 

* Ibid. , p.  21. 
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a correlation  should  exist  between  a coherence  estimate  near  unity  and 
an  accurate  value  of  the  calculated  directional  angle.  This  correla- 
tion is  the  basis  for  using  the  coherence  parameter  to  evaluate  the 
accuracy  of  the  calculated  directional  angles  at  each  frequency. 


Sources  of  Error  in  Calculating  Target  Direction 


21.  Application  of  the  directional  algorithm.  Equation  11,  is 
susceptible  to  several  sources  of  error.  Generally,  these  errors  can 
be  separated  into  two  categories:  (a)  instrumentation  and  signal- 

processing errors,  and  (b)  wave  propagation  errors.  Errors  in  the 
first  category  arise  from  phase  response  characteristics  of  the  instru- 
mentation, finite  spatial  size  of  the  geophones  (i.e.,  they  do  not 
sample  a point  in  space  as  assumed  in  the  derivation  of  the  directional 
algorithm) , and  numerical  errors  of  a statistical  nature  that  occur 
within  the  signature  processing  techniques.  Errors  in  the  second  cate- 
gory arise  from  refraction  and  reflections  of  the  seismic  wave  on  its 
path  from  the  target  to  the  geophones  and  from  errors  in  the  plane 
wave  approximation  (paragraph  11) . These  errors  are  further  treated 
in  Appendix  A.  An  additional  potential  source  of  error  arises  from 
the  difficulty  in  proper  interpretation  of  the  times  and  f^c 

from  the  calculated  phase  shifts  AQ^  and  A0^,  (Equations  6 and  7). 
The  coherent-signature  analysis  techniques  (paragraph  20)  mathematically 
constrain  the  magnitude  of  the  calculated  phase  shifts  between  -it 
and  it  , even  though  the  actual  phase  shift  may  be  the  calculated  value 
+2nr  or  -?.utt,  where  n is  an  integer.  This  mathematical  ambiguity  of 
the  true  phase  shift  is  an  expression  of  a basic  physical  difficulty. 

In  Figure  6,  two  sinusoidal  time  series  are  presented  that  are  shifted 
in  phase  relative  to  each  other.  The  question  is:  Which  peaks  in  one 

series  are  associated  with  which  peaks  in  the  other?  The  ambiguity 
problem  can  be  eased  by  shortening  the  distance  between  the  geophones 
in  the  array.  However,  an  ambiguity  will  always  occur  when  the  magni- 
tude of  the  phase  shift  exceeds  it  . Since  the  phase  shift  is  propor- 
tional to  the  frequency  (Equation  6) , there  exists  a cutoff  frequency 
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above  which  the  ambiguity  problem  always  occurs.  The  ambiguity  can  be 
removed  for  frequencies  greater  than  the  cutoff  frequency  by  adding  or 
subtracting  2nn  to  the  magnitude  of  the  calculated  phase  shift,  where  n 
is  the  appropriate  integer. 


Figure  6.  Comparison  of  two  time  series  shifted  in  phase 

relative  to  each 
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PART  III:  APPLICATION  OF  SIGNAL-PROCESSING 

TECHNIQUES  TO  SIGNATURES  OF  MOVING  TARGETS 


Data  Collection 


22.  Two  test  programs  were  conducted  to  obtain  data  for  evaluating 
the  seismic  target-location  system.  The  instrumentation  and  geophone- 
array  deployment  were  tested  in  two  geologically  dissimilar  test  sites, 
and  the  results  were  subjected  to  analysis  by  the  techniques  described 
in  paragraphs  19  and  20. 

23.  The  first  test  program  was  conducted  in  January  1974  at 

Vicksburg,  Mississippi.  The  second  was  conducted  in  April  1974  at  Fort 
Carson,  Colorado.  The  tests  were  designed  around  the  following  vari- 
ables: vehicle  type  (i.e.  M35,  M.113,  M151,  and  M728) , vehicle  velocity, 

vehicle  number  (i.e.  one  or  tv/o  vehicles  moving  simultaneously),  surface 
condition  (i.e.  on-road  or  cross-country),  range,  geological  conditions 
(i.e.  Vicksburg  and  Fort  Carson  test  sites),  array  spacing  (2  m and 

6 m) , and  geophone  type  (Scientific  RS-10  and  Mark  Products  L-4  geo- 
phones). The  configuration  of  each  of  the  short-based  arrays  used  in 
these  two  test  programs  is  presented  in  Figure  7.  Test  variables  for 
both  test  programs  are  tabulated  in  Appendix  C. 

Vickburg  tests 

24.  The  Vicksburg  test  site  consisted  of  deep  lean  clay  (loess) 
soils  with  a water  table  about  30  cm  below  the  surface.  Vehicles 
tested  were  an  M113  APC  (light  tracked  vehicle),  an  M35  (2-1/2-ton 
truck),  and  an  M151  (jeep).  The  layout  of  the  Vicksburg  test  site  is 
shown  in  Figure  8.  The  geophone  array,  denoted  by  the  triangular 
symbol,  consisted  of  six  vertical  Geospace  HS-10  geophones.  Each 
geophone  was  implanted  deep  enough  to  allow  its  upper  surface  to  be 
covered  by  several  centimetres  of  top  soil  to  reduce  the  possibility  of 
acoustic  signals  coupling  with  the  geophones.  The  data  from  each 
geophone  were  simultaneously  recorded  on  a 14-track  tape  recorder  at 
tape  speeds  of  3-3/4  and  7-1/2  ips.  At  either  tape  speed,  the  frequency 
response  of  the  recorder  is  well  above  1000  Hz.  The  HS-10  geophone  is 
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Figure  8.  Vicksburg  test  site  layout 
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Road  B 


sensitive  to  frequencies  from  1 Hz  to  somewhat  greater  than  350  Hz.  In 
the  single-target  tests,  the  vehicles  were  run  on  roads  A and  B (Figure 
8),  beginning  at  the  indicated  starting  points  and  proceeding  to  the 
intersection  of  the  two  roads.  Each  test  was  terminated  slightly  before 
the  vehicle  reached  the  intersection  because  a constant  velocity  could 
not  be  maintained  while  negotiating  the  turn.  The  position  of  the 
vehicle  was  monitored  by  placing  event  markers  at  50-m  intervals  along 
the  respective  roads.  Radio  contact  was  maintained  during  each  test  run 
so  that  a passenger  in  the  vehicle  could  identify  the  location.  By 
recording  the  time  at  which  each  event  marker  was  passed,  and  inter- 
polating to  locate  the  position  between  event  markers,  an  accurate 
record  of  the  vehicle  location  was  obtained  throughout  the  test. 

Passage  of  event  markers  was  recorded  for  each  vehicle  over  radio  links 
of  two  different  frequencies.  A number  of  trial  runs  were  made  prior  to 
the  actual  test  sequence  so  that  drivers,  passengers,  and  instrument 
operators  could  become  familiar  with  the  recording  procedure.  It  was 
found  that  the  vehicle  position  at  the  stakes  could  be  easily  monitored 
to  within  5 m or  less  using  the  radio  communication  link. 

Fort  Carson  tests 

25.  The  geology  of  the  Fort  Carson  test  site  was  complex.  It  con- 
sisted of  a 20-  to  30-cm  layer  of  silt  overlying  fractured  limestone. 
Vehicles  tested  were  the  M113  APC,  the  M35,  the  M151,  and  an  M728  (heavy 
tracked  combat  engineering  vehicle  mounted  on  an  M60  chassis).  The 
layout  of  the  test  site  is  shown  in  Figure  9.  The  vehicles  were  run 
across  country  (path  A)  and  on  a compacted  surface  (path  B,  a well -used 
tank  trail).  The  Fort  Carson  tests  were  conducted  in  essentially  the 
same  manner  as  in  the  Vicksburg  tests.  Vehicles  were  run  on  two  inter- 
secting paths,  A and  B,  with  straight-line  runs  of  500  and  550  m,  respec- 
tively, as  indicated  in  Figure  9.  Two  arrays  of  geophones  were  deployed 
in  the  manner  indicated  in  Figure  7 with  a spacing  between  the  arrays  of 
220  m.  The  most  westward  array  consisted  of  triaxial  Mark  Product  L-4 
geophones;  the  other,  of  Scientific.  HS-10  geophones.  Data  were  recorded 
simultaneously  for  both  arrays. 
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Results  of  Signature  Analysis 


26.  The  two  techniques  of  signature  analysis,  cross-correlation 
and  coherent-signature  (paragraphs  19  and  20) , were  applied  to  data 
selected  from  the  two  test  programs.  Only  those  data  (see  Appendix  C) 
required  to  demonstrate  the  applicability  of  the  algorithms  were  ana- 
lyzed because  of  time  and  fund  constraints.  No  rigorous  criteria  were 
used  to  select  data  from  specific  tests,  nor  were  data  selected  so  as  to 
make  either  procedure  give  better  results,  i.e.,  the  data  selected  were 
typical  of  those  recorded  at  each  test  site  and  were  not  anomalous  in 
any  way.  The  results  of  the  application  of  these  techniques  to  the 
selected  data  are  discussed  below. 

Results  of  cross-correlation  analysis 

27.  Cross-correlation  analysis  was  applied  to  data  selected  from 
the  Vicksburg  tests  as  listed  below.  Test  numbers  correspond  to  those 
in  Table  Cl. 


Test 

No. 

Vehicle(s) 

Velocity 

km/hr 

Closest 

Range 

m 

Farthest 
Range,  m 

Road 

(Figure  8) 

Array 

Spacing,  m 
(Figure  7) 

1 

Ml  13 

32 

130 

600 

B 

2 

2 

Mil  3 

32 

130 

600 

B 

6 

4 

M35 

24 

70 

280 

A 

6 

8 

M3  5 

32 

130 

410 

B 

6 

M151 

24 

78 

210 

A 

These  tests  were  selected  for  analysis  because  they  include  both  wheeled 
and  tracked  vehicles  traveling  at  representative  velocities  and  usefully 
distant  ranges,  in  addition  to  a multiple-target  test.  Directional 
angles  were  calculated  at  various  ranges  between  the  range  bounds 
identified  in  the  above  tabulation.  The  vehicles  were  " tracked"  along 
their  respective  roads  (A  and  B in  Figure  8).  Data  were  analyzed  for 
the  M113  tests  using  both  the  2-  and  6-m  geophone  arrays  identified  in 
Figure  7.  In  each  of  the  tests,  the  seismic  signatures  were  digitized 
at  a rate  of  1500  samples/sec.  The  continuous  integral  of  Equation  12 
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reduces  to  a discrete  summation  over  the  range  of  integration  (i.e.  -T/2 
to  T/2).  The  increment  in  the  time  domain,  dt  , is  then  given  by 
(1/1500)  sec.  An  integration  range  defined  by  T ■ 0.2  sec  was  found 
to  be  acceptable  for  analysis  because  the  period  of  the  seismic  waves  in 
the  dominant  energy  band  was  less  than  that  value.  The  results  of  the 
data  analysis  are  discussed  below. 

28.  Test  1.  Test  1 was  conducted  with  a 2-m  geophone  array  and  an 
M113  APC  moving  at  32  km/hr.  Measured  and  calculated  directional 
angles  as  the  M113  moved  down  road  B are  presented  in  Table  1;  the 
average  absolute  error  is  8.3  deg.  The  measured  and  calculated  direc- 
tional angles  are  plotted  as  a function  of  vehicle  travel  time  in 
Figure  10a.  Time  was  chosen  as  a parameter  rather  than  range  in  this 
figure,  because  range  does  not  uniquely  determine  the  position  of  the 
vehicle  in  the  test  geometry.  The  position  of  the  vehicle,  although  not 
tracked  continuously,  is  described  by  a continuous  line.  The  corre- 
sponding range  to  the  vehicle  as  a function  of  time  can  be  obtained  from 
Figure  10b. 

29.  Test  2.  Test  2 was  conducted  with  a 6-m  geophone  array  and  an 
M113  APC  moving  at  32  km/hr.  Measured  and  calculated  directional  angles 
as  the  M113  moved  down  road  B are  presented  in  Table  1;  the  average 
absolute  error  is  4.4  deg.  The  measured  and  calculated  directional 
angles  are  plotted  as  a function  of  vehicle  travel  time  in  Figure  11a. 
The  corresponding  range  to  the  vehicle  as  a function  of  time  can  be 
obtained  from  Figure  lib. 

30.  Test  4.  Test  4 was  conducted  with  a 6-m  array  and  an  M35  2- 
1/ 2-ton  truck  moving  at  24  km/hr.  Measured  and  calculated  directional 
angles  as  the  M35  moved  down  road  A are  presented  in  Table  2;  the 
average  absolute  error  is  4.7  deg.  The  measured  and  calculated  direc- 
tional angles  are  plotted  as  a function  of  vehicle  travel  time  in 
Figure  12a.  The  corresponding  range  to  the  vehicle  as  a function  of 
time  can  be  obtained  from  Figure  12b. 

31.  Test  8.  Test  8 was  conducted  with  a 6-m  array  and  two  vehi- 
cles, an  M35  truck  and  an  M151  jeep  moving  at  32  and  24  km/hr,  respec- 
tively. Measured  and  calculated  directional  angles  as  the  M35  moved 


31 


RANGE,  M OtRECTtONAL  ANCLE,  OEG 


b.  MEASURED  RANGE  FROM  SENSOR  ARRAY  TO  TARGET 

Figure  10.  Position  information  for  an  M113  APC  (velocity 
- 32  km/hr,  geophone  spacing  -2m)  using  geophones  k , 5, 

and  6 (Figure  J) 
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RANGE,  M DIRECTIONAL  ANGLE-  DEG 


d.  CALCULATED  AND  MEASURED  DIRECTIONAL  ANGLES 


b.  MEASURED  RANGE  FROM  SENSOR  ARRAY  TO  TARGET 

Figure  11,  Position  information  for  an  M113  APC  (velocity 
= 32  km/ hr,  geophone  spacing  = 6 m)  using  geophones  1,  2, 

and  3 (Figure  7 ) 
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PANGE,  M DIRECTIONAL  ANGLE.  DEG 


dL.  CALCULATED  AND  MEASURED  DIRECTIONAL  ANGLES 


b.  MEASURED  RANGE  FROM  SENSOR  ARRAY  TO  TARGET 


Figure  12.  Position  information  for  an  M3  5 truck  (velocity 
= 2b  km/hr,  geophone  spacing  = 6 m)  using  geophones  1,  2, 

and  3 (Figure  7) 
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down  road  B and  the  M151  moved  down  road  A are  presented  in  Table  3.  As 
pointed  out  in  paragraph  16,  the  geophones  sense  the  superposition  of 
the  two  sets  of  signatures  emanating  from  the  two  targets.  Whenever  the 
dominant  bands  (i.e.  bands  containing  the  bulk  of  the  energy  from  the 
target)  overlay  each  other  so  that  the  individual  signatures  of  the  two 
targets  compete  for  the  calculated  directional  angle,  the  analysis  of 
the  signatures  results  in  a calculated  directional  angle  that  is  likely 
to  be  somewhere  between  the  two  directions  to  the  vehicles.  Comparison 
of  measured  and  calculated  directions  in  Figure  13  indicates  that  the 
calculated  directional  angles  lie  generally  between  the  two  measured 
directions  to  the  vehicles.  However,  the  calculated  directional  angles 
are  closer  to  the  direction  of  the  M35  than  of  the  M151,  probably 
because  the  signature  from  the  former  vehicle  generally  dominated  that 
emanating  from  the  lighter  M151  vehicle.  The  absolute  errors  were 
calculated  relative  to  the  measured  location  of  the  M35;  the  average  was 
20.3  deg. 

32.  Summary . Based  on  the  cross-correlation  analysis  of  the 
selected  data,  the  following  comments  can  be  made  concerning  the  use  of 
short-based  arrays  for  calculating  directional  angles i 

a.,  The  array  with  6-m  spacing  can  be  used  to  obtain  direc- 
tional angles  within  an  average  accuracy  of  less  than  5 deg 
in  a single-target  (or  quiet-background)  environment.  This 
comment  should  be  tempered  by  the  fact  that  the  accuracy  of 
the  method  is  terrain  dependent,  and  these  results  are 
derived  from  only  one  test  site. 

_b.  The  accuracy  of  the  technique  decreases  considerably  in 
a multiple-target  environment  relative  to  a single-target 
environment. 

£.  The  accuracy  of  the  technique  decreases  for  the  2-m  array 
relative  to  the  6-m  array.  Reasons  for  this  decrease  in 
accuracy  as  the  geophone  spacing  is  shortened  are  discussed 
in  Appendix  A. 

The  limitations  of  cross-correlation  analysis,  particularly  the  require- 
ment of  a large  spacing  for  useful  accuracy  (i.e.  6-m)  and  the  decrease 
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RANGE,  M DIRECTIONAL  ANGLE.  DEG 


a.  CALCULATED  AND  MEASURED  DIRECTIONAL  ANGLES 


b.  MEASURED  RANGE  FROM  SENSOR  ARRAY  TO  TARGET 

Figure  13.  Position  information  for  multiple-target  test  in- 
volving M35  truck  and  Ml 51  Jeep  (velocity  = 2h  km/hr,  geo- 
phone spacing  = 6m)  using  geophones  1,  2,  and  3 (Figure  7) 
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in  accuracy  in  a mil  1 1 T ple-target  environment,  suggest  that  alternative 
techniques  cf  signature  analysis  should  be  investigated.  Coherent- 
signature  analysis  was  investigated  in  order  to  address  the  shortcomings 
of  the  simpler  cross-correlation  analysis  (paragraph  17). 

Results  of  coherent-signature  analysis 

33.  The  coherent-signature  analysis  was  applied  to  data  selected 
from  the  Vicksburg  and  Fort  Carson  tests  as  listed  below.  Test  numbers 
correspond  with  those  in  Table  Cl. 


Test 

No. 

Site 

Vehicle (s) 

Velocity 

km/hr 

Range , m 

Road 

or 

Path 

Array 

Spacing 

m 

Geophone 

Type 

1 

Vicksburg 

Ml  13 

32 

450 

B 

2 

HS-10 

5 

Vicksburg 

M3  5 

16 

185 

A 

2 

HS-10 

9 

Ft.  Carson 

M728 

16 

280 

A 

2 

L-4 

7 

Vicksburg 

M35 

32 

190 

B 

2 

HS-10 

M151 

24 

100 

A 

These  data  were  selected  because  they  include  signatures  from  wheeled 
(M151,  M35) , light  tracked  (M113) , and  heavy  tracked  (M728)  vehicles; 
they  include  representative  ranges  and  velocities  and  data  from  geologi- 
cally different  test  sites;  and  two  geopnone  types  (i.e.  vertical-sensing 
Scientific  HS-10  and  triaxial-sensing  Mark  Product  L-4)  are  represented. 
(Only  the  vertical  sensing  component  of  the  triaxial  geophones  was  used 
in  the  analysis.)  The  2-m  spacing  between  geophones  in  the  array  was 
used  in  <..11  tests  instead  of  the  6-m  spacing  because  the  smaller  array 
is  less  susceptible  to  the  ambiguities  in  the  computed  phase  shift 
(paragraph  21),  particularly  at  the  higher  frequencies.  As  in  the 
cross-correlation  analysis,  the  data  were  digitized  at  1500  samples/sec, 
allowing  a frequency  analysis  up  to  750  Hz.  The  results  of  the  d ita 
analysis  are  discussed  below. 

34.  Test  1.  Test  1 was  conducted  at  the  Vicksburg  site  with  an 
Mll^  moving  at  32  km/hr  on  road  B at  a range  of  450  m.  The  array 
consisted  of  HS-10  geophones  spaced  2 m apart.  Calculated  directional 
angles  for  the  M113  at  a measured  directional  angle  of  68.5  deg  are 
presented  in  Table  4.  The  corresponding  frequencies  and  coherence 
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estimates  are  also  given  in  the  table.  For  low  frequencies  (less  than 
10  Hz)  the  calculated  directional  angles  are  consistently  too  small 
in  spite  of  the  high  coherence  values.  This  discrepancy  in  the  cal- 
culated directional  angles  was  also  typical  of  the  results  of  analysis 
of  the  remaining  tests  and  could  be  due  to  different  phase  response 
characteristics  of  the  geophones  at  the  low  frequencies.  Above  10  Hz, 
the  accuracy  of  the  calculated  directional  argles  is  closely  related  to 
the  coherence  value.  The  calculated  angles  are  consistently  accurate 
for  coherence  values  above  0.85,  although  some  calculated  angles  are 
quite  accurate  (i.e.  within  5 deg)  for  low  coherence  values,  such  as 
the  angle  at  114.37  Hz.  The  average  of  the  calculated  directional 
angles  having  coherence  values  greater  than  0.85  and  frequency  above 
10  Hz  is  69.5  deg,  which  compares  very  favorably  with  the  measured  value 
of  68.5  deg.  A polar  plot  showing  the  calculated  directional  angles 
for  each  frequency  greater  than  10  Hz  and  having  coherence  estimates 
greater  than  0.85  is  presented  in  Figure  14. 

35.  Test  5.  Test  5 was  conducted  at  the  Vicksburg  test  site  with 
an  M35  moving  at  16  km/hr  on  road  A at  a range  of  185  m.  The  array 
consisted  of  HS-10  geophones  spaced  2 m apart.  Calculated  directional 
angles  for  the  M35  at  a measured  directional  angle  of  248  deg  are 
presented  in  Table  5 for  those  frequencies  greater  than  10  Hz  and  having 
associated  coherence  estimates  greater  than  0.85.  A polar  plot  of  the 
calculated  angles  appearing  in  Table  5 is  presented  in  Figure  15. 

36*  Test  9.  Test  9 was  conducted  at  the  Fort  Carson  site  with  an 
M728  moving  at  16  km/hr  on  path  A (cross-country)  at  a range  of  280  m. 
The  array  consisted  of  Mark  Product  L-4  geophones  spaced  2 m apart. 
Calculated  directional  angles  for  the  M728  at  a measured  directional 
angle  of  311  deg  are  presented  in  Table  6 for  those  frequencies  greater 
than  10  Hz  and  having  associated  coherence  estimates  greater  than  0.85. 

A polar  plot  of  the  calculated  directional  angles  appearing  in  Table  6 
is  presented  in  Figure  16. 

37.  Test  7.  Test  7 was  conducted  at  the  Vicksburg  test  site  with 
two  vehicles,  an  M35  and  an  M151,  moving  at  32  and  24  km/hr  on  roads  B 
and  A,  respectively.  The  array  consisted  of  HS-10  geophones  spaced  2 m 
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particular  frequency  for 
coherency  >0.85.  Frequency 
is  proportional  to  distance 
from  inner  circle  (0  Hz). 
Outer  circle  is  200  Hz. 


irection  to  an  M113  at 
t site. 
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0# 


(array  position)  Outer  circle  is  200  Hz. 


Figure  15.  Polar  plot  of  direction  to  an  M35  at 

Vicksburg  test  site 
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Measured  position 
for  M728  at  1 6 km/hr 
and  range  of  280  m / 


Scale  of  frequency  - relative  to  center 

(array  position) 


Figure  16.  Polar  plot  of 
Fort  Carson 
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Note:  Dot a represent  calculated 

directional  angle  at  a 
particular  frequency  for 
coherency  >0.85.  Frequency 
is?  proportional  to  distance 
from  inner  circle  (0  Hz). 
Outer  circle  is  200  Hz. 


apart.  Calculated  directional  angles  for  location  of  the  M35  at  a 
measured  angle  of  45  deg  and  the  M151  at  a measured  angle  of  236  deg  are 
presented  in  Table  7 for  those  frequencies  greater  than  10  Hz  and  having 
an  associated  coherence  estimate  greater  than  0.85.  A polar  plot  of  the 
calculated  directional  angles  appearing  in  Table  7 is  presented  in 
Figure  17.  Three  clusters  of  calculated  directional  angles  appear.  The 
first  cluster  consists  of  four  points  whose  average  directional  angle  is 
238.4  deg.  This  average  calculated  directional  angle  compares  favorably 
with  the  measured  angle  for  the  M151  jeep  (236  deg).  The  second  cluster 
consists  of  two  calculated  angles  at  approximately  17  deg.  No  target 
was  known  to  be  in  that  direction.  The  anomaly  could  be  due  to  coherent 
instrument  noise,  coherent  background  noise  (i.e.  another  seismic  source 
or  target),  or  the  random  nature  of  the  coherence  estimate  (i.e.,  the 
process  of  estimating  coherence  from  random  data  is  itself  random  because 
estimates  based  upon  finite  time  intervals  have  a probability  of  being 
in  error) . The  third  cluster  consists  of  a series  of  nine  calculated 
angles  between  the  frequency  range  of  115  to  176  Hz.  The  average 
directional  angle  for  this  cluster  is  45.5  deg,  which  agrees  favorably 
with  the  measured  angle  for  the  M35  (45  deg) . 

38.  Summary.  Based  upon  the  coherent-signature  analysis  of  the 
selected  data,  the  following  comments  can  be  made  concerning  the  use  of 
the  short-based  array  for  calculating  directional  angles: 

£.  The  short-based  array  of  2-m  spacing  can  be  used  to  obtain 
directional  angles  well  within  an  average  absolute  error  of 
5 deg  in  a single-target  (quiet-background)  environment. 

This  statement  is  supported  by  results  obtained  at  two 
geologically  dissimilar  sites  (i.e.  Vicksburg  and  Fort 
Carson) . 

b^.  The  coherent-signature  analysis  appears  to  be  a potentially 
useful  approach  for  tracking  multiple  targets.  Although 
two  spurious  estimates  occurred  in  analysis  of  the  multiple- 
target  data  (paragraph  37),  in  an  operable  system  such 
"targets"  could  be  studied  over  sequential  time  intervals 
to  determine  if  they  are  generated  by  randomness  in  the 
directional  angles,  by  coherent  background  noise  in  the 
array  environment,  or  by  a "target"  having  characteristics 
(i.e.  velocity,  frequency,  range,  etc.)  of  the  targets  of 
interest. 
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c.  Seismic  signatures  produced  by  wheeled  and  tracked  vehicles 
contain  sufficient  energy  between  0 and  200  Hz  to  permit 
computation  of  directional  angles  at  ranges  at  least  as 
large  as  450  m (Test  1,  paragraph  34) . 


N 

0° 


Figure  17.  Polar  plot  of  directions  to  multiple  targets 
consisting  of  an  M35  and  an  M151 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


39.  Based  upon  the  study  reported  herein,  the  following  conclusions 
can  be  drawn: 

a.  A short-based  array  of  geophones  can  be  used  to  obtain 

““  data  for  determining  usefully  accurate  directional  angles 
to  single,  moving  targets  by  cross-correlation  analysis 
and  coherent-signature  analysis  (paragraphs  32  and  38, 
respectively) . Coherent-signature  analysis  does  not 
require  that  the  geophone  spacing  be  as  large  as  that 
required  for  equivalent  accuracy  using  cross-correlation 
analysis  (2  m versus  6 m for  cross-correlation  analysis 
using  the  HS-10  geophones) . It  must  be  recognized  that 
the  limits  on  accuracy  as  the  spacing  in  the  array  is 
changed  depend  upon  the  spatial  dimensions  of  the  geophones 
(i.e.  18  cm  for  a HS-10)  and,  therefore,  smaller  arrays 
could  probably  be  developed,  provided  smaller  transducers 
are  placed  in  the  array. 

b.  Coherent-signature  analysis  is  a potential  technique  for 
tracking  multiple  targets  (paragraph  37).  Cross- 
correlation analysis  does  not  result  in  acceptably  accurate 
values  of  the  directional  angle  in  a multiple-target 
environment  (paragraph  31) . 

_c . The  seismic  location  system  reported  herein,  together  with 
its  associated  signal-processing  techniques  used  to 
calculate  directional  angles,  does  not  appear  to  be 
exceedingly  sensitive  to  terrain  variations.  The  extent 
and  limitations  of  terrain  variations  require  further 
investigation. 

Recommendations 


40.  It  is  recommended  that: 

a.  Investigation  of  the  seismic  techniques  of  target  location 
be  continued.  The  data  collected  during  the  field  program 
described  herein  (Vicksburg  and  Fort  Carson  programs)  should 
be  analyzed  to  better  determine  the  sensitivity  of  target 
location  methods  to  terrain  variation.  Furthermore,  as 
the  data  on  hand  are  not  extensive,  additional  data  should 
be  collected  at  geologically  dissimilar  sites,  and  these 
data  should  be  analyzed  as  described  above. 
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b_.  Instrumentation  response  characteristics,  including  those 
of  small  seismic  transducers,  be  investigated  to  determine 
the  state-of-the-art  capability  of  constructing  a short- 
based  array  of  minimum  spatial  dimensions.  The  sensitivity 
of  the  directional  algorithm  to  changes  in  design  criteria 
could  be  evaluated  as  discussed  in  Appendix  A. 

The  study  be  continued  to  determine  the  feasibility  of 
incorporating  target  classification  and  location  into  the 
capabilities  of  a single  geophone  array.  The  utility  of 
combined  acoustic  and  seismic  sensors  for  ranging 
(Appendix  B)  and  classification  should  be  included;  in 
particular,  the  feasibility  of  applying  the  ranging 
technique  (paragraph  7b)  to  firing  artillery  should  be 
further  investigated.  Alternative  techniques  of  signal 
processing,  in  addition  to  those  reported  herein,  that 
make  use  of  phase  shift  and  amplitude  properties  of  the 
signatures  should  be  considered. 

<1.  A study  be  initiated  to  define  the  performance  envelope 
for  competing  or  supplementing  techniques,  such  as 
acoustic  techniqvies , as  a function  of  terrain  and 
vegetation  cover  in  order  to  compare  (or  supplement)  seismic 
and  acoustic  approaches  to  target  location  efficiently. 
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Table  1 

Position  Information  for  M113  APC  Using 
2-  and  6-m  Arrays 


Directional  Angles,  deg 

Absolute 

Time,  sec 

Range,  m 

Calculated 

Measured 

Error,  deg 

16.6 

410 

2-m  Array 
10 

15 

5 

22.5 

365 

11 

17 

6 

28.5 

315 

19 

19 

0 

34.4 

270 

24 

22 

2 

40.3 

215 

45 

25 

20 

45.9 

190 

36 

29 

6 

51.7 

155 

42 

36 

6 

57.5 

130 

13 

44 

29 

63.5 

130 

59 

52 

7 

70.0 

145 

72 

74 

2 

Avg  8.3 

0 

600 

6-m  Array 
10 

12 

2 

5.23 

510 

12 

13 

1 

16.6 

410 

14 

15 

1 

22.5 

365 

14 

17 

3 

28.5 

315 

17 

19 

2 

34.4 

270 

31 

22 

9 

40.3 

215 

39 

25 

14 

45.9 

190 

34 

29 

5 

51.7 

155 

41 

36 

5 

57.5 

130 

48 

44 

4 

63.5 

130 

58 

52 

6 

70.0 

145 

73 

74 

1 

Avg  4.4 

Table  2 

Position  Information  for  M3 5 Using  6-m  Array 


Directional  Angle,  deg 

Absolute 

Time,  sec 

Ranee,  m 

Calculated 

Measured 

Error,  dej 

15.9 

280 

245 

256 

11 

29.5 

230 

256 

254 

2 

44.4 

180 

243 

250 

7 

60.2 

135 

241 

242 

1 

74.4 

90 

231 

228 

3 

90.1 

70 

200 

197 

3 

105.3 

105 

159 

154 

5 

120.4 

145 

134 

128 

6 

Avg  4.7 

Table  3 

Position  Information  for  Multiple  Targets  Using  6-m  Array 


M3  5 

Location 

M151 

. Location 

Calculated 

Time 

sec 

Range 

m 

Directional 
Angle,  deg 

Range 
m 

Directional 
Angle,  deg 

Directional 
Angle,  deg 

Absolute 
Error.* **  dec 

18.3 

410 

15 

210 

253 

6 

9 

25.3 

3b5 

17 

165 

248 

59 

42 

31.3 

315 

19 

131 

241 

36 

17 

38.1 

270 

22 

93 

230 

55 

33 

43.5 

215 

25 

78 

212 

32 

7 

48.9 

190 

29 

80 

185 

44 

15 

55.3 

155 

36 

107 

150 

60 

24 

62.4 

130 

44 

144 

138 

56 

12 

68.3 

130 

52 

a* 

** 

70 

18 

74.2 

145 

74 

** 

** 

58 

26 

Avg  20.3 

* Error  is  computed  using  M35  directional  angles  because  cross- 
correlation analysis  technique  results  in  computed  angles  that 
’'attempt"  to  track  only  the  M35. 

**  M151  completed  passage  on  road  A (Figure  8)  at  62,7  sec  into  test 

run.  Vehicle  continued  moving  for  several  seconds  after 
completion  time. 


Table  4 

Directional  Angles  Calculated  for  M113  APC  at 
Measured  Angle  of  68.5  deg,  Using  2-m  Array 


Frequency,  Hz 

Coherence 

Calculated  Directional 
Angle,  deg 

2.93 

0.97 

63.2 

4.39 

0.99 

51.4 

5.86 

0.99 

52.2 

7.33 

0.99 

52.8 

8.79 

0.99 

53.3 

10.26 

0.99 

63.9 

11.73 

0.97 

70.0 

13.20 

0.95 

74.5 

14.66 

0.92 

69.6 

16.13 

0.82 

60.2 

17.59 

0.68 

52.6 

19.06 

0.73 

64.8 

20.52 

0.68 

49.9 

21.99 

0,61 

25.4 

23.46 

0.61 

22.2 

24.92 

0.62 

134.5 

26.39 

• 

0.60 

• 

73.5 

• 

• 

• 

68.92 

• 

• 

0.46 

• 

• 

86.2 

70.38 

0.51 

252.6 

71.85 

0.21 

275.1 

73.31 

0.40 

5.8 

74.78 

0.43 

11.4 

76.25 

0.81 

66.9 

77.71 

0.91 

73.5 

(Continued) 


Table  4 (Concluded) 


Frequency,  Hz 


Coherence 


Calculated  Directional 
Angle,  deg 


107.04 

0.50 

289.2 

108.51 

0.57 

290.2 

109.97 

0.79 

295.3 

111.14 

0.90 

66.3 

112.91 

0.90 

68.3 

114.37 

0.65 

72.9 

115.84 

0.46 

78.0 

Table  5 

Directional  Angles  Calculated  for  M3 5 at 


Measured  Angle  of  248  deg,  Using  2-m  Array 


Frequency , Hz 

Coherence 

Calculated  Directional 
Angle,  deg 

10.26 

0.98 

249.6 

11.73 

0.99 

249.9 

13.19 

0.98 

250.2 

14.66 

0.98 

247.5 

16.13 

0.99 

247.2 

19.06 

0.97 

246.5 

20.53 

0.91 

245.3 

21.99 

0.92 

248.6 

23.46 

0.95 

250.7 

24.92 

0.96 

250.7 

26.39 

0.97 

248.7 

27.86 

0.97 

246.1 

29.32 

• 

0.97 

• 

245.7 

• 

• 

• 

158.36 

• 

• 

0.86 

• 

• 

247.6 

173.02 

0.91 

247.4 

174.49 

0.93 

247.3 

178.89 

0.91 

248.4 

180.36 

0.92 

247.5 

181.82 

0.89 

245.5 

Table  6 

Directional  Angles  Calculated  for  M728  at 


11.72 

0.99 

315.8 

13.67 

0.98 

313.8 

15.62 

0.95 

316.2 

76.17 

0.89 

305.4 

80.01 

0.86 

314.2 

82.03 

0.91 

314.3 

83.98 

0,93 

314.1 

85.94 

0.92 

313.7 

97.66 

0.85 

304.3 

128.91 

0.85 

310.8 

130.86 

0.88 

310.4 

132.81 

0.90 

309.5 

134.77 

0.86 

308.4 

140.63 

0.87 

306.7 

142.58 

0.90 

307.2 

144.53 

0.97 

308.5 

146.48 

0.98 

308.5 

148.44 

0.98 

308.1 

150.39 

0.98 

307.6 

152.34 

0.97 

307.3 

156.25 

0.94 

308.1 

158.20 

0 1 88 

308.3 

160.16 

0.92 

309.2 

162.11 

0.95 

309.6 

164.06 

0.^5 

309.0 

166.02 

0.95 

308.0 

(Continued) 


Table  6 (Concluded) 


Frequency t Hz 


Coherence 


Calculated 

Directional  Angle,  de% 


167.97 

0.87 

307.3 

193.36 

0.92 

310.8 

195.31 

0.95 

309.6 

197.27 

0.92 

309.6 

199.22 

0.92 

308.4 

Table  7 

Directional  Angles  Calculated  for 
Multiple-Target  Tests,  Using  2-m  Arra 


Calculated 


Frequency,  Hz 

Coherence 

Directional  Angle,  deg 

26.39 

0.88 

236.2 

27.86 

0.88 

241.4 

32.36 

0.89 

236.0 

33.72 

0.89 

240.0 

86.51 

0.96 

16.6 

87.99 

0.95  * 

18.0 

115.84 

0,90 

43.9 

117.30 

0.90 

42.9 

130.50 

0.90 

44.8 

131.97 

0.92 

44.8 

145.16 

0.87 

44.8 

146.63 

0.89 

45.8 

173.02 

0.95 

48.0 

174.50 

0.97 

47.6 

175.96 

0.94 

46.6 

APPENDIX  A:  SENSITIVITY  OF  THE  DIRECTIONAL  ALGORITHM 

TO  TEMPORAL  AND  SPATIAL  ERRORS 

1.  The  directional  algorithm  for  using  a short-based  array  of 

gecphones  is  susceptible  to  several  sources  of  error.  Generally,  these 
errors  can  be  separated  into  two  categories:  (a)  instrumentation  and 

signal-processing  errors,  and  (b)  wave  propagation  errors.  Errors  in 
the  first  category  are  usually  well  defined  and  include  errors  arising 
from  the  spatial  response  of  the  geophones  (i.e.,  transducers  sample  an 
area  of  space  rather  than  a point).  Errors  in  the  second  category  are 
generally  defined  by  empirical  analysis.  However,  an  analysis  of  prop- 
agation errors  due  to  deviations  from  the  plane  wave  approximation 
(i.e. , the  wave  is  more  accurately  described  in  two  dimensions  by  a 
circle)  can  be  analytically  examined. 

Instrumentation  and  Signal-Processing  Errors 

* 

2.  The  algorithm  for  using  an  array  of  geophones  placed  at  the 
vertices  of  an  equilateral  triangle  is  given  by  Equation  11  of  the  main 
text  as 


3.  Assuming  either  that  the  phase  response  of  the  instrumenta- 
tion is  specified  by  a known  value  or  that  the  accuracy  of  the  signal- 
processing technique  for  determining  the  phase  shift  between  the  signa- 
tures is  defined,  the  sensitivity  of  the  algorithm  to  errors  in  xAT, 

AB 

and  tAC  can  determined.  If  At  bounds  the  uncertainty  in  the 
measurement  of  x^,  and  x^  , the  resulting  uncertainty  A<j>  in  the 
directional  angle  <J>  can  be  obtained  by  differentiating  Equation  Al. 

The  resulting  expression  for  A<f>  is  bounded  for  reasonably  small  values 
of  Ax  by 


Al 


A<J>  < — -r—*  F(4>)  * (57.29  deg/radian)  (A2) 

with 

F (4>)  “ (VT  + tan  <j>)  cos  for  0 c|>  <_  60°  (A3) 

where  v is  the  velocity  of  the  seismic  wave  at  a given  frequency,  & 
is  the  geophone  spacing,  and  A<j>  is  expressed  in  degrees.  The  angular 
dependence  of  the  uncertainty,  F(<j>)  , is  graphically  presented  in 
Figure  A1  and  falls  within  the  range  1.73  to  2.  The  angular  dependence 
of  F (q> ) in  Equation  A3  is  defined  only  in  the  region  0 <_  <J>  <_  60°. 
However,  by  symmetry  arguments,  this  function  can  be  repeated  through 
the  entire  range  of  <f)  . 

4.  As  shown  in  Figure  A1  the  uncertainty  A<J>  , which  is  pro- 
portional to  F(<f>)  , is  smallest  when  the  incoming  wave  approaches  along 
one  of  the  lines  formed  by  any  two  geophones  in  the  array  (i.e.  along 

AC  or  AB)  and  the  largest  when  the  wave  approaches  from  an  angle  bisecting 
these  lines.  It  should  be  recognized  that  this  analysis  will  result  in  a 
worst-case  description  of  the  uncertainty  A<}>  because  the  errors  in 
and  are  superimposed  in  the  most  destructive  manner. 

5.  The  inequality  (Equation  A2)  can  be  used  to  obtain  design 

parameters  in  the  following  manner.  Suppose  that  the  design  criteria 
are  that  A<|>  be  less  than  5 deg  and  the  array  dimensions  be  less  than 
0.5  m.  Suppose  further  that  the  wave  velocity  in  the  operational 
environment  is  400  m/sec  in  the  frequency  range  of  interest.  Since 
F(<j>)  <2,  the  inequality  can  be  used  to  define  the  accuracy  requirement 
on  the  time  measurements  (i.e.  and  t^,)  as  being  less  than  50 

microseconds  (50  x 10  ^ sec).  The  response  time  of  the  instrumentation 
system  could  then  be  specified  as  being  less  than  50  microseconds  for 
frequencies  up  to  400  Hz.  The  cutoff  frequency  of  400  Hz  is  convenient 
because  ambiguities  in  the  calculated  phase  do  not  occur  below  that 
value  (i.e.,  the  spacing  in  the  array  is  less  than  half  the  wavelength 
for  frequencies  below  400  Hz) . 


A2 


Figure  Al.  Angular  dependence  of  the  uncertainty  in  the  directional 
measurement,  F(<j>).  F(<|>)  is  defined  by  Equation  A3. 


A3 


Errors  Due  to  Plane  Wave  Assumption 


6.  As  described  previously,  the  derivation  of  Equations  10  and 
11  in  the  main  text  depends  on  the  assumption  that  the  approaching  waves 
are  plane  waves.  They  are  actually  described,  in  two  dimensions,  by 
circular  wave  fronts,  if  the  medium  is  assumed  reasonably  isotropic.  The 
magnitude  of  the  resulting  error  arising  from  the  plane  wave  assumption 
can  be  investigated  by  considering  the  size  of  the  first  correction 
term  arising  from  linear  dependence  of  the  x^g  and  x^c  upon  the 
ratio  of  the  array  spacing  to  the  range  to  the  target  l/x  . To  this 
order  of  approximation  (i.e.  linear  dependence  of  x^g  and  x^c  upon 
i/r),  the  following  expression  can  be  derived  for  the  equilateral  array. 


TAC  1 . /J  . , , i sin2  <J>  - sin2  (60°  - <J>) 

— - Y + r tan  * *27  L cos  ♦ 


(M) 


7.  The  last  term  of  Equation  A4  is  the  first-order  correction 
term  to  the  ratio  of  Xw,  and  x^g  . The  sensitivity  of  the  algorithm 
to  this  correction  can  be  examined  by  differentiating  Equation  A1  with 
respect  to  the  ratio  tac/tab  * resulting  uncertainty,  A<J>  , in 

degrees  is  then  bounded  by 

A<|>  < (0.933  ~)  x (57.29  deg/radian) 

Therefore,  for  ranges  greater  than  11  *:imes  the  array  spacing,  the 
uncertainty  in  <J>  is  less  than  5 deg.  In  most  practical  cases,  errors 
resulting  from  the  use  of  the  plane  wave  assumption  are  insignificant 
for  short-based  arrays  (i.e.  less  than  7-m  spacing). 


A4 


APPENDIX  B:  TARGET  LOCATION  USING  THE  DIFFERENCES  IN 

ARRIVAL  TIMES  OF  SEISMIC  AND  ACOUSTIC  WAVES 


1.  Reliance  on  triangulation  techniques  for  target  location  has 
serious  drawbacks  in  a battlefield  environment,  because  they  require 
that  at  least  two  groups  of  geophones  be  emplaced.  For  this  reason,  it 
may  be  desirable  to  investigate  ways  by  which  both  direction  and  range 
to  the  target  could  be  obtained  from  one  geophone  array.  For  example, 
one  possible  approach  for  determining  the  range  to  impulsive  energy 
sources,  such  as  artillery  firings  or  impacts,  is  to  measure  the 
differences  in  arrival  times  of  seismic  and  acoustic  waves.  Acoustic 
signatures  are  usually  obtained  with  microphones;  however,  geophones 
near  the  ground  surface  are  sensitive  to  both  seismic  and  acoustic 
energy.  This  approach  is  limited  to  impulsive  sources,  because 
correlation  between  acoustic  and  seismic  energy  for  continuous  sources 
is,  at  best,  extremely  difficult.  The  concept  of  ranging  in  this  manner 
depends  on  differences  between  the  velocities  of  the  acoustic  and  seismic 
waves.  A geophone  signature  resulting  from  the  firing  of  a 155~mm 
howitzer  is  presented  in  Figure  Bl.  Initially,  the  signature  amplitude 
is  small,  arising  primarily  from  random  background  noise.  Then  the 
signature  amplitude  increases,  subsides,  and  finally  rises  sharply  to 
saturation.  A simultaneous  acoustic  recording  at  the  location  of  the 
geophone  revealed  that  the  saturated  portion  of  the  geophone  signature 
occurred  in  the  same  time  interval  as  the  passage  of  the  acoustic  wave 
from  the  firing  of  the  howitzer.  Therefore,  the  saturated  portion  of 
the  geophone  signature  is  attributed  to  the  acoustic  energy.  The  lesser 
energy  region  preceding  the  acoustic  energy  is  attributed  to  the  seismic 
energy  emanating  from  the  firing.  The  leading  edge  of  this  seisraically 
active  region  is  not  well  defined.  However,  it  precedes  the  time  of 
arrival  of  the  acoustic  signature  by  approximately  0.5  sec.  A technique 
needs  to  be  developed  to  accurately  identify  the  arrival  time  of  the 
seismic  signature.  Such  a technique  could  be  based  upon  cross  correla- 
tion of  the  signatures  received  at  the  array. 
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2.  A geophone  array  can  be  used  to  determine  the  direction  to 
the  target  by  application  of  the  directional  algorithm  (Equation  11  of 
the  main  text).  The  velocity  of  the  seismic  wave  train  can  be  obtained 
from  the  expression 


v(u>) 


&AC  cos  <t> 


T 


AC 


which  is  a rearrangement  of  Equation  9 of  the  main  text.  The  range  to 
the  target,  r , can  then  be  obtained  from  the  expression 


r 


sc 


V 7 

s a 


v _ 

s 


— At 
a 


where 

vq  * velocity  of  the  seismic  wave 

vq  “ velocity  of  the  acoustic  wave 

At  “ difference  in  time  between  the  initial  receptions  of  the 
acoustic  and  seismic  waves 

The  velocity  of  the  acoustic  wave  in  air  is  accurately  approximated  by 
330  m/sec  in  windless  environments.  For  maximum  accuracy  the  wind 
velocity  would  have  to  be  included  in  the  value  of  v . 

3.  The  approach  described  above  is  potentially  capable  of 
providing  a method  of  placing  counterbattery  fire  on  a target  weapon, 
particularly  if  that  weapon  is  in  a fixed  location.  The  procedure  is 
hypothesized  as  follows: 

a.  The  location  of  the  enemy  weapon  is  estimated  by  cal- 
culating the  direction  of  arrival  of  its  seismic  signa- 
ture, using  the  three-geophone  array  previously  described. 
The  range  is  obtained  from  the  difference  between  the 
arrival  times  of  the  acoustic  and  seismic  waves,  as 
described  in  the  previous  paragraph  (Figure  B2) . 

b.  A counterbattery  weapon  is  then  fired  at  the  estimated 
location  (Figure  B3) , and  the  location  of  the  impact  of 
that  round  is  calculated,  using  the  geophone  array.  Since 
the  propagation  paths  will,  in  most  cases,  be  very  similar, 
it  can  be  assumed  that  the  errors  in  the  directional  angle 
and  range  will  be  essentially  the  same  as  the  errors  in 
the  calculated  location  of  the  enemy  weapon  relative  to 
its  actual  location. 
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True  direction 
to  enemy  weapon 


Error  in  directional  estimate 


Direction  calculated  using 
geophone  array 


LEGEND 


True  location  of  enemy  weapon 


Initial  estimated  location  of 
enemy  weapon 


Figure  B2.  Initial  target-location  analysis 
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CEP  of  counter battery  weapon 


Figure  B3.  Calculation  of  position  of  first-round  impact 

of  counterbattery  fire 
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£.  The  location  of  the  counterbattery  impact,  as  calculated 
by  the  geophone  array,  is  compared  with  the  location  of 
the  target,  as  calculated  in  the  same  way.  Since  the 
actual  location  of  the  impact  within  the  central  error  of 
probability  (CEP)  is  unknown,  it  is  assumed  that  the  impact 
of  the  counterbattery,  round  is  exactly  on  the  calculated 
target  position.  The  difference  between  the  first  estimate 
of  the  location  of  the  impact  of  the  counterbattery 
weapon  and  the  assumed  location  of  the  first  impact  is 
taken  as  an  estimate  of  the  magnitude  of  the  range  and 
azimuth  correction  required  for  the  counterbattery  weapon 
to  hit  the  target.  The  range  and  azimuth  corrections 
required  to  presumably  hit  the  enemy  weapon  are  in  the 
opposite  direction  of  the  errors  themselves.  In  Figure  B3 

Range  correction  ■ - r^ 

Azimuth  correction  “ $2  “ 4>i 

d..  With  these  corrections,  a second  counterbattery  round  is 
fired,  as  illustrated  in  Figure  B4.  The  estimated  location 
of  the  second  round,  as  calculated  from  the  data  from  the 
geophone  array,  can  now  be  compared  with  the  initial 
estimate  of  the  hostile  weapon  location.  A new  cycle  of 
range  and  azimuth  corrections  would  still  further  refine 
the  position,  so  that  subsequent  rounds  could  be  fired  for 
effect. 

4.  Although  the  signature  analysis  techniques  presented  in  this 
report  have  been  applied  only  to  moving  targets  (i.e.  APC's,  trucks, 
etc.),  th^re  is  no  barrier  to  extending  one  of  these  (cross-correlation 
analysis)  to  firing  artillery  and  impact  signatures.  The  directional 
algorithm.  Equation  11  of  the  main  text,  has  been  applied  to  firing 
artillery  and  impacts,  using  graphic  techniques  that  simulate  cross- 
correlation analysis.*  However,  the  coherent-signature  analysis  tech- 
nique, as  presented  in  this  study,  is  not  applicable  to  impulsive 
targets,  such  as  firing  artillery,  because  of  the  requirement  for  a time 
average  in  the  analysis  scheme. 


* WESFE  letter  dated  27  Aug  73  to  Project  Manager,  REMBASS , subject: 
"Seismic  Methods  for  Target  Location." 
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Figure  BU . Calculation  of  range  and  azimuth  corrections 
(from  Figure  B3)  for  second  round  of  counterbattery  fire 
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APPENDIX  C:  TABULATION  OF  SINGLE-  AND  MULTIPLE-TARGET  TEST  CONDITIONS 


1.  Table  Cl  lists  the  test  conditions  for  single-  and  multiple- 
target  data  collected  in  this  program.  Ranges  from  the  geophone  array 
varied  continually  between  the  approximate  limits  of  100  m and  600  m. 
(Exact  ranges  can  be  identifed  from  the  data  logs  and  the  geometry  of 
the  event  markers  relative  to  the  geophone  array.  See  paragraphs  23,  24, 
and  25  of  the  main  text  for  a discussion  of  the  test  procedures.)  Geo- 
phone arrays  consisted  of  the  geophone  types  identified  in  the  tabulation. 
Spacing  of  geophones  in  the  arrays  is  described  in  Figure  7 of  the  main 
text. 


Cl 


Vehicle  Test  Data 


(Continued)  

* Hard  surface  is  a well-used,  hard  packed  road.  Natural  surface  is  cross-countr r travel. 

**  Second  target  is  the  associated  multiple  target.  NA  indicates  that  the  test  w as  a single-target 
test. 


Table  Cl  (Concluded) 


In  accordance  vith  BR  70-2-3#  paragraph  6c(l)(b), 
dated  15  February  1973#  a facsimile  catalog  card 
in  Library  of  Congroes  format  is  reproduced  belov. 
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